Piezoresistive heat engine and refrigerator 
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Heat engines provide most of our mechanical power and are essential for transportation on macro- 
scopic scale. However, although significant progress has been made in the miniaturization of elec- 
trostatic engines, it has proven difficult to reduce the size of liquid or gas driven heat engines below 
10 7 /xm 3 . 

Here we demonstrate that a crystalline silicon structure operates as a cyclic piezoresistive heat 
engine when it is driven by a sufficiently high DC current. A 0.34 fim 3 engine beam draws heat from 
the DC current using the piezoresistive effect and converts it into mechanical work by expansion 
and contraction at different temperatures. This mechanical power drives a silicon resonator of 
1.1 xlO 3 fxia into sustained oscillation. Even below the oscillation threshold the engine beam 
continues to amplify the resonator's Brownian motion. When its thermodynamic cycle is inverted, 
the structure is shown to reduce these thermal fluctuations, therefore operating as a refrigerator. 



Most heat engines operate by the cyclic expansion and 
contraction of a gas at different temperatures. Much ef- 
fort has gone into miniaturizing these conventional heat 
engine concepts* — . However this downscaling is diffi- 
cult, because the efficiency and power density of these 
engines reduce strongly with its size*. A promising route 
for the miniaturization of heat engines is the use of an 
electrically heated solid as a working substance, because 
the large heat capacitance of solids allows much higher- 
power densities and electrical interconnect facilitates en- 
ergy transport and localized heat generation at the mi- 
croscale. This concept has often been applied in elec- 
trothermal actuators for microscopic devices 5 - - — , how- 
ever the similarity between thermal expansion actuators 
and heat engines has received little attention. 

In order to operate as a cyclic heat engine, a thermal 
actuator needs to be supplemented with a mechanism 
that regulates the thermodynamic cycle of heating, ex- 
pansion, cooling and compression. This cycle can be reg- 
ulated by an external AC frequency synthesizer— or by 
a feedback loop, where a sustained reciprocating motion 
is generated from a DC powe r 5 ' 9 ' 10 . In previous reports 
of electrothermal engines, the thermodynamic cycle was 
always regulated by external electronic transistor circuits 
and amplifiers. 

Here it is shown that a well-dimensioned silicon crystal 
of uniform composition operates as a cyclic heat engine 
when driven by a DC current. It is demonstrated that 
a silicon engine beam can drive a silicon resonator into 
sustained oscillation. In contrast to earlier work 5 -!***, no 
external transistors or amplifiers are needed to regulate 
the heat engine, because a mechanism based on the in- 
trinsic piezoresistive heating and thermal expansion of 
silicon provides the required feedback. The engine dis- 
placement volume of the heat engine, which similar to car 
engines is defined by the space occupied by the working 
substance, is 0.34 /im 3 . Despite its low efficiency, the en- 
gine's power density is almost a factor 1000 higher than 
that of modern car engines. 

Like any other reciprocating heat engine, this piezore- 
sistive heat engine requires a minimum threshold power 
Pt to compensate its internal dissipation and sustain its 



cyclic motion. In the first part of this article the opera- 
tion of the piezoresistive heat engine above this threshold 
power P t is demonstrated by driving a resonator into sus- 
tained oscillation. In the second part the motion of the 
device at driving powers below P t is investigated. Mea- 
surements show that the engine beam continues to in- 
crease the thermal energy stored in the resonance mode 
like a heat pump. When the thermodynamic cycle of 
the engine is reversed, it is shown to reduce the stored 
thermal energy in the resonant mode, similar to recently 
developed optical and RF refrigerators^—. In contrast 
to conventional refrigerators, which simultaneously re- 
move thermal energy from many degrees of freedom of 
a substance, this type of refrigerators cools by removing 
thermal energy from only a single mechanical resonance 
mode. This piezoresistive refrigerator provides an alter- 
native way to reduce undesirable Brownian fluctuations 
in micromechanical sensors and mirrors. 



Piezoresistive heat engine 

Experiments are performed on the homogeneous crys- 
talline silicon structure shown in Fig. [1^,. The struc- 
ture consists of a mass measuring 12.5x60.0x1.5 /im 3 , 
which is suspended by a 3 fira wide spring beam and a 
280 nm narrow engine beam. Both beams have a length 
of 800 nm. Before discussing the stand-alone operation 
of the heat engine, its in-plane mechanical resonance at 
1.26 MHz (Fig. [5^,) is characterized by actuating it via 
terminal T3 with an AC electrostatic force and detect- 
ing the strain variations in the engine beam using the 
piezoresistive effect (sec Methods). To operate the de- 
vice as a piezoresistive heat engine, the external AC ac- 
tuation voltage is disconnected from terminal T3. Thus, 
as shown in Fig. [TJd the device is only connected at ter- 
minal Tl to a DC current source Ida and to a capacitively 
coupled oscilloscope. All other terminals are grounded. 
At low values of Ida only noise is observed on the oscil- 
loscope. However, if the DC current is increased above 
a threshold 7 f =1.19 mA a remarkable effect occurs: the 
device spontaneously starts to oscillate and generates a 
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FIG. 1: Micrograph, schematics (not to scale), measurement circuit and operation mechanism of the piezore- 
sistive heat engine and refrigerator, a Scanning electron microscope (SEM) image of the heat engine. The inset shows a 
magnification of the wide spring beam and narrow engine beam by which the resonator mass m is suspended, b Schematic of 
the device and electrical connections used for the measurements in Fig. [5p,c,d. c Schematic of the device and electrical circuit 
used for the measurements in Fig. [3] and [4] d Besides the elastic coupling between displacement and force, a thermally delayed 
feedback mechanism via the internal electrical and thermal variables occurs in a piezoresistive spring in the presence of a DC 
current flow, e Schematic thermodynamic cycle in the piezoresistive spring, when the device is operating at a single frequency 
as a heat engine. As a result of thermal delay, the temperature Ti, and stress a are phase shifted with respect to the strain e 
and piezoresistive heating power p pr . The axes indicate the longitudinal strain e and stress cr, which are proportional to the 
volume change V and pressure P respectively. Therefore the engine beam expands at a higher temperature than at which it is 
compressed. It thus generates an amount of mechanical work W e from piezoresistive heat during each cycle, f Illustration of 
the phases CD-© of the thermodynamic cycle shown in which are identified by the AC resistive heating power p pr and beam 
temperature Ti,. To illustrate the position dependent heating power a fictitious external heat sink (blue) and source (red) are 
drawn. In reality the internal heating power in the beam depends on its position as a result of the piezoresistive effect. 



sinusoidal output voltage v ac with a frequency of 1.26 
MHz. Oscilloscope and spectrum analyser measurements 
of the signal are shown in Fig. [^Jd, and[5Jl. The motion 
of the device is also studied by stroboscopic illumination 
under a microscope. The stroboscope is triggered by the 
electrical output v ac of the heat engine. By adjusting the 
phase delay of the stroboscope a slow-motion movie of 
the engine's sinusoidal motion is made (see Supplemen- 
tary Video). Stroboscopic and electrical measurements 
show that just above It the displacement amplitude of 
the mass depends very sensitively on the DC current Idc- 
If Idc is increased more than 50 /iA above It , the optically 
observed mechanical amplitude x stops increasing and is 
saturated by collisions of the resonator mass against elec- 
trode T3. This amplitude saturation is also observed in 
the piezoresistance variations r ac . 

The observed oscillation of the heat engine is at- 
tributed to the thermodynamic feedback mechanism 
which is schematically shown in Fig. QJl. When an al- 
ternating mechanical displacement x = x$e lu}t with fre- 
quency uj is applied to a piezoresistive beam with spring 



constant fc, this will not just result in an elastic force, 
but will also result in a resistance variation r ac via the 
piezoresistive effect. If a constant DC current Idc is flow- 
ing through the spring, the resistive heating power there- 
fore changes by p pr = l\ c r a c and causes a temperature 
change, which in turn will generate a thermal expansion 
force F te that adds to the external forces F ext on the 
spring, such that kx = F te + F ext . 

The amplitude and the thermal phase-delay of the 
thermal expansion feedback force F te can be described by 
a complex coefficient /3 which is defined by F te = fH\ c kx. 
The coefficient /3 can be calculated from the material con- 
stants and geometry of the spring (sec Supplementary 
Discussion A). The response of the piezoresistive spring 
to an external force can now be described by a complex 
spring constant k* s = F ext /x = k(l — (31^)- When the 
piezoresistive spring is coupled to a mass m, it forms a 
resonator that can store energy E/j similar to the flywheel 
of a conventional heat engine. In the absence of current 
(Idc = A), the resonator's intrinsic damping can be ex- 
pressed by a Q-factor Q; nt and its resonance frequency 
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FIG. 2: Heat engine output signal, a Characterization of the in-plane mechanical bending resonance by external excitation 
at different chamber pressures. b,c,d Operation of the heat engine at Id c =h+0. 01 mA at a chamber pressure P=0.01 mbar. 
In panel b the spontaneous start-up of the engine is measured by an analogue oscilloscope when a current 7d c =1.20 mA is 
switched on at t—0 s. The self-sustained sinusoidal output signal is measured c by a digital oscilloscope and d by a spectrum 
analyser that determines the power spectral density (PSD). From the oscilloscope data in panel c, the amplitude xo of the 
centre-of-mass is estimated to be a;o=43 nm, using xq = v ac .o/(VdcK pr ), with Vd c =0.99 V and an AC amplitude iw,o=27 mV. 



is given by oj = ^Jk/m. In the small-signal approx- 
imation, the linearised equation of motion of the heat 
engine driven resonator is equivalent to the well-known 
harmonic oscillator: mx + + k* s x = 0. For a si- 

nusoidal motion x = xoe lult near the resonance frequency 
ojq the imaginary (damping) part of this equation be- 
comes zero at a threshold current value I?=(Q i nt Im f3)~ 1 . 
Above this DC current Idc > It the power generated by 
the heat engine beam exceeds the intrinsic losses and the 
resonator will be brought into sustained oscillation by the 
engine beam. 

The heat capacitance of the beam results in a thermal 
delay between piezoresistivc heating power and temper- 
ature in the engine beam. As a consequence, the beam 
expands at a higher beam temperature Xj, than at which 
it contracts. The amplitude of temperature variations 
during the oscillation measured in Fig. [5J: is estimated 
to be 0.2 K by finite element simulations. During each 
period of oscillation the silicon working substance in the 
piezoresistive beam therefore goes through a thermody- 
namic cycle, similar to a Stirling engine, during which 



it converts part of the resistive heat into work W e . The 
piezoresistive beam therefore operates as a heat engine. 

Figure []J shows the P — V diagram and Fig. QJ the 
phases of the thermodynamic cycle of the engine. In 
phase ® the speed of the mass results in an expansion 
of the engine beam at high temperature. In phase © 
the strain in the beam is tensile, such that it is cooling 
because the piezoresistance of n-type silicon is low. As 
a consequence the temperature of the beam is reduced 
in phase ® such that the beam contracts at low tem- 
perature, until it reaches phase @ where the compressive 
strain causes a high piezoresistive heating power, result- 
ing again in an elevated temperature in phase ©. 



Heat pump and refrigerator 

After having shown the operation of the heat engine 
above the oscillation threshold It , we now investigate its 
characteristics at lower DC currents. Below the threshold 
the displacement x of the resonator will still exhibit ran- 
dom thermally excited 'Brownian' motion. The spectral 
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FIG. 3: Voltage spectral density <v^> for different values of Idc- Measurements by a spectrum analyser using the 
circuit in Fig. []J show that a in the open state (Zi oa d = oo) the thermal motion is amplified when the heat engine is driven 
by a current Idc- b For Z\ OBiC \ — 40/JH, the mechanical energy decreases with increasing current and the device operates as a 
refrigerator. Note the factor 4 difference in frequency scale, c 1/T e g as determined from the data in[3£i and [3b by evaluating 
^fcsTcff = \k <x 2 >, where <x 2 > is determined by integration of the fitted dashed lines (see Methods). The difference between 
the j/-axis intersection of the data and the ambient temperature is within experimental uncertainty limits. 



displacement noise density <x 2 J > of this motion near the 
resonance frequency in the presence of the feedback force 
F te can be determined from the fluctuation-dissipation 
theore m 11 ' 14 ! 17 : 



<x 2 „> 



2^k B T b /{kQ 

int ) 

{u 2 - ul) 2 + (wtJo/Qoff) 2 

Vint 



(1) 

(2) 



Where k B is Boltzmann's constant and it is assumed that 
Qeff > 1 and \(3I 2 C \ < 1. 

To determine <x 2 J >, the AC voltage spectrum <v 2 >, 
which is proportional to <x 2 J >, is measured (Fig. |3^,b) 
at several values of Idc by a spectrum analyser using 
the measurement circuit in Fig. QJ. In Fig. [3]d an in- 
ductor Zi oa d = 40/iH is placed in parallel to the struc- 
ture. Equation (JTJ) is used to fit the data (see Meth- 
ods). The measurements in Fig. [3ji show that even be- 
low the threshold current It the thermal fluctuations in 
the resonance mode are amplified. It is observed that 
the stored thermal energy in the resonance mode in- 
creases with DC current, such that it equals the energy 
which would be stored in the mode if the resonator would 
be at an effective temperature T e g which is defined by 
\k B T cK = \k f°° <x 2 > ^ = hk <x 2 >. Similar to res- 

onators cooled by optical or RF methods^— the tem- 
perature T e g represents the mechanical energy stored in 
the resonance mode corresponding to the degree of free- 
dom x and is not representative for the other degrees 
of freedom of the resonator. By integrating the fitted 
<x 2 J > curves, T e g is determined and is plotted in Fig. 
|3fc. Depending on the load impedance Zi oa d, the effective 



temperature T e g of the resonance mode either increases 
like in a heat pump (Fig. |3fi) or is refrigerated down 
to an effective temperature of 70 K at /dc=3.0 mA (Fig. 
[3b). 

The difference between the spectra in Figs. [3^ and 
(5}d is caused by the fact that the feedback coefficient j3 
is proportional to the ratio = P P r/{Id c r ac) between 
the piezoresistive heating power p pr and the AC resis- 
tance r ac as can be seen from Fig. QJL Therefore the 
coefficient ft can be written as the product of two com- 
plex numbers: /3 = 7zA)j where /?o depends solely on 
the device geometry and material parameters and jz 
can be controlled externally via Z\ oac \ in the following 
way. Figure [TJ; shows that a resistance change r ac of the 
beam will induce both an AC current i ac — ~v ac /Z to t 
and AC voltage v ac = i?dcW + Idcfac- The total heat- 
ing power Pdc + Ppr in the engine beam is given by 



(Idc+lacY(R be 

that 7z = 1 — 2R be 



r ac ) and for r ac <C Rbe 

m/ (Ztot + Rdc)- 



it follows 



Quantitative analysis 

Voltage spectra like those in Fig. |3Ji,b are recorded for 
several real and imaginary values of Z\ oac i, using resistors, 
capacitors and inductors. The fitted values of Q e g are 
plotted in Fig. @Jt and b. For small values of Zi oa d a more 
accurate determination of Q e g was made by generating 
a random electrostatic noise force on the resonator using 
a white voltage noise source connected to terminal T3. 
The observed linear dependence of 1/Q ft on ij c in Fig. 
3Ji,b corresponds well with equation @. From the slope 
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FIG. 4: Measurements and fits to confirm the operation mechanism proposed in Fig. [T]i. a,b Inverse effective 
Q-factor 1 /Q e ff versus DC current l\ c for different values of the load impedances Z\ oss) \ indicated in the graph. Dashed lines are 
fits using equation (fJJ). Symbols on the a;-axis correspond to the minimum current level at which self-sustained oscillation was 
observed, a Resistive load impedances Z\ 0& &. b Capacitive and inductive load impedances Zi oa ,d- c Measured slopes — Im /3 of 
the curves in Fig. a,b plotted along the real (left) and the imaginary (right) 7z-axis for resistive (red circle), capacitive (blue 
square) and inductive (green diamond) values of Zi oa d. An excellent fit of the data is obtained by plotting — Im (yzPo) with 
po=— 123.7 + 64. 6iA -2 , for resistive (dashed lines), capacitive and inductive (solid lines) values of Z\ ^. If Im /3 is positive the 
engine beam acts as a heat pump that increases the effective temperature and if Im /3 is negative the engine beam acts as a 
refrigerator that cools the effective temperature of the resonance mode. 



of these linear fits, Im j3 is determined and is plotted 
against 72 as symbols in Fig. Hp. As shown by solid 
and dashed lines in Fig. @Ji, an excellent multiple linear 
regression fit of the data is obtained using the function 
Im ft = Im ("fz A)), which yields the fit parameters (3q = 
-123.7 + 64.6i A' 2 and R beam = 439.3 ft. 

In Supplementary Discussion A an analytical model 
is derived for the complex Young's modulus and spring 
constant k* s and it is used to derive an estimate for /?o, 
which yields /3q = —132 + 157i A~ 2 . A finite element 
simulation of the full geometry including the anisotropy 
of the silicon crystal results in /3q = — lll+58z A~ 2 . Both 
the excellent fits in Fig. [3] and |U and the quantitative 
agreement between the measured and simulated values 
of /3 support the proposed feedback mechanism in Fig. 

m. 



Flow rates of heat and work 

In order to gain insight in the operation of the heat en- 
gine, a schematic showing the rate equations for the flow 
of heat and work between electrical source, heat sink, en- 
gine beam and resonance mode is shown in Fig. [5^,. All 
rates shown in this figure are time averaged, such that all 
variations which cancel each other during a single cycle 
are averaged out. The three modes of operation which 
were experimentally observed in Fig. [5fi, [5fx, and [3p are 
shown in Fig. [Sp. For Im /3 < 0, Whe is negative, thus 
cooling the effective temperature of the resonance mode 
and operating as a refrigerator. For < Qint^rf c Im/? < 1 
heat is pumped into the resonance mode, increasing T c ff 
until a new steady state is reached at which W e + Ui = 0. 
In this mode the engine beam operates as a heat pump. 
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FIG. 5: Diagrams showing the heat flow and operation modes of the engine beam, a Rate equations governing 
the energy transport between heat sink, engine beam and resonance mode. The resistive heating Q r powers the engine beam, 
which converts the heat into work Whe on the resonance mode. The resonance mode dissipates the power Whe at a rate 
Qfd and uses the excess power to perform external work W e or to increase T e g. Heat from the engine beam flows away to 
the heat sink by its thermal conductance ku- Internal energy (Uib,Ui) is stored by the heat capacitances of the beam at 
constant volume (C„b) and of the resonance mode (|fcs). All quantities in this diagram are time averaged values and the dots 
represent time derivatives. Integration of equation |T]) yield a 11 ' 14 ' 17 the relation T e s/Tb = Qcff/Qint and combining this with 
equation Q and the definition of Qi nt results in the equations for Qfd and Whe- b Three different regimes of operation are 
distinguished. Depending on the feedback factor Qi nt /| c Im/?, the device operates as a refrigerator, heat pump or heat engine. 
If no external work or heat is transferred (We + Ui=0), the effective temperature of the resonance mode follows the blue line 

Toff/Tft = (1 - Qint/ilm/3)" 1 = Qeff/Qint. 



For Qi nt I% c Imf3 > 1, the steady state rate equation in 
Fig. [5^i can only be satisfied with W e + Ui > 0. There- 
fore external work W e is generated and the engine beam 
operates as a heat engine. The distinction between the 
work and heat flows in Fig. [5] is subtle. For T e ff — > oo it 
follows that Qoff tends to infinity such that the motion of 
the resonator is perfectly sinusoidal, like in Fig. This 
energy can theoretically be stored with 100% efficiency 
and should be called work. At any finite value of T e g the 
resonator performs a thermally driven random motion, 
which cannot be fully converted into work according to 
Carnot's theorem^ and should be called heat. 



Performance 

For the miniaturization of engines a high power den- 
sity Pdens is essential, since it determines the minimum 
size the engine can have in order to perform its task 
within the specified time. The maximum power Whe that 
the demonstrated heat engine can generate is given by 
ifc < x 2 > cjo^ c Im/3. In order to achieve the maximum 
power density, a current Idc=5 mA is applied for which 
the amplitude of the centre of mass is limited by colli- 
sions with electrode T3 to x=100 nm. For an optimum 
value of Im /3=100 A -2 (sec Fig. |4j:) it is found that 
W} le =l2.7 nW, about 4% of which is dissipated by the 
intrinsic losses (Qfd)- It thus follows that the piezoresis- 
tive heat engine, with an engine displacement volume of 
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the engine beam of V = 0.34 /zm 3 , has a power density 
Pbcam = W e /V=37 GW/m 3 which is almost a factor 1000 
higher than that of car engines which have typical power 
densities of around 50 MW/m 3 . The power density Pbeam 
of the piezoresistive heat engine is also high when com- 
pared to other microengine o 1 ' 20 ' 21 , which have reported 
power densities up to 3 GW/m 3 . However, in some of 
these cases the power density is defined using the vol- 
ume of the total structure instead of the engine displace- 
ment volume. Due to the relatively large size of the proof 
mass compared to the engine beam size, this definition 
of power density is less favourable for the piezoresistive 
heat engine, yielding p tot =ll MW/m 3 . 



Using the equations from Fig. [5ji, the efficiency of 
the presented engine is calculated to be r]he = (Whe ~ 
Qfd)/Qr =6.4xl0 -7 . Besides being limited by the 
Carnot theorem, the efficiency of the piezoresistive heat 
engine is also limited by the fracture limit of silicon to 
around 7? pr ,si,max = 3.4 x 10 -4 , as analysed in the Supple- 
mentary Discussion B. It is therefore expected that the 
efficiency and power density of the piezoresistive heat en- 
gine can still be improved by a factor ^500 by increasing 
the displacement amplitude and thermal design of the 
engine. Even after optimization, this efficiency is low 
compared to that of micro- and nano-engincs driven by 
electrostatic forces^ 2 -"— . Nevertheless, the comparison 
in Supplementary Discussion B shows that piezoresistive 
heat engines have the potential to outperform electro- 
static engines in terms of power density. This high power 
density and the absence of transistors can facilitate the 
development of mechanical actuator— and senso r 10 ! 27 ar- 
rays with very high element densities. 



Similar to other recently investigated cooling 
methods^—, the coefficient of performance of the 
refrigerator 77 coo i = —Whe/Qr is very small when T e g 
is at or below room temperature. The refrigerator will 
therefore not be useful for lowering the temperature 
of systems with many degrees-of-freedom. However, 
because the heat capacity -§jr^ of a single-degree-of- 

freedom resonance mode is only \k,B and because the 
thermal insulation of the mode is high MS (X result of its 
high Qmt (see Fig. [5£i), the piezoresistive refrigerator 
can still reduce T e g by more than 200 K. The Brownian 
fluctuations associated with the energy i/csT e ff often 
dominate the noise in mechanical sensors and mirrors 
and their reduction can thus increase the signal-to-noisc 
ratio and might ultimately enable quantum-limited 
measurement and control of mechanical resonators. 
The presented piezoresistive refrigerator provides an 
alternative cooling method for reducing these Brownian 
fluctuations, which has the advantage that it does not 
require laser or RF sources and only requires a silicon 
chip with a relatively low DC power. 



Methods 

Structure The piezoresistive heat engine in Fig. [T]i 
is made out of a 1.5 /im thick silicon layer on a silicon- 
on-insulator wafer, with a phosphor doping concentra- 
tion A^=4.5xl0 18 cm" 3 giving a specific resistivity of 
/9 f / c =10" 4 f2m. This doping concentration results in a 
low temperature coefficient of resistivity- 2 ^ of less than 
0.1% per K (see Supplementary Fig. C). The thin crys- 
talline silicon layer is structured in a single mask step 
by a deep reactive ion etch and the buried Si02 layer 
below the mass and beams is removed in a hydrogen flu- 
oride vapour etch. A DC current Ida is driven through 
the beams via terminals Tl and T2 as shown in Fig. [TJd 
and c. As a result of the geometry of the structure, the 
current density, heating power density and mechanical 
stress are concentrated in the narrow engine beam, but 
the resonance frequency of the resonator is mainly de- 
termined by the spring beam and resonator mass. The 
in-plane mechanical bending resonance mode that de- 
termines the fundamental resonance frequency /o=1.26 
MHz has a spring constant fc=256 N/m as determined 
from finite element method (FEM) simulations. Mea- 
surements are performed at 315 K in a vacuum chamber 
at a pressure below 10 -2 mbar. 

External excitation The in-plane fundamental me- 
chanical resonance of the structure at 1.26 MHz is char- 
acterized by actuating it via terminal T3 with an AC 
electrostatic force F ac ^ act = eo^Kic.actWac^ct/s 2 gener- 
ated by a voltage Vdc.act + v aCt3 _ ct on terminal T3, with 
Vdc,act= -1 V. The actuation gap c/=200 nm has an ac- 
tuation area A = 60 x 1.5 /im 2 . A DC current /^c = 0.1 
mA is applied as in Fig. [Ud. The displacement of the 
centre of mass x is proportional to the measured piezore- 
sistive AC current i pr , which is shown in Fig. [5^. The 
reduction of the intrinsic Q-factor Qj nt upon increasing 
the pressure in the vacuum chamber confirms the me- 
chanical nature of the resonance. The solid fits in Fig. 
[2k allow us to determine the piezoresistive factor K pr , 
which is defined by r ac / Rdc = K pr x, where r ac is the 
amplitude of the piezoresistance variation and Rdc is the 
total DC resistance. From the maximum piezoresistive 
current i pr ^ pea k the piezoresistive factor is determined to 
be = k = _g 6 x 1Q 5 m -i 

Spectrum analysis The spectra in Fig. [3] and 2] are 

measured using the electrical circuit in Fig. [TJ:. Be- 
cause the displacement of the centre-of-mass x u = x(ui) 
causes a piezoresistive voltage v u = v ac = [Z^~ l + 
R^]~ 1 Id c K pr x (Fig. [Tb), the electrical spectrum <i' 2 > 
measured by the spectrum analyser is proportional to 
equation JTJ). The voltage spectral density in Fig. [3^ 
and b is fit by dashed curves which are given by equation 
© with <vl>= lj c K^<xl>/[Z^ + i?j/] 2 and which 
include a fit parameter to account for the white John- 
son noise background. The total AC impedance parallel 
to the resonator Z tot = (iwCo)" 1 + (iujCi + Z^^) -1 
can be controlled by adjusting ^i oa d- The DC resis- 
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tance R^c = -Rbeam + i?o=794 ft, bias and cable capac- 
itances Co = 54 nF and C\ = 70 pF were measured 
using an impedance analyser at 1.26 MHz. Terminals 
T3 and T4 are grounded. A 1 MO input impedance 
voltage buffer amplifier is used to minimize the effects 
of cables and spectrum analyser on the engine opera- 
tion. In Fig. [3Jd DC resistive heating can lead to a tem- 
perature increase of the beam which opposes the refrig- 
eration mechanism. It is estimated from the tempera- 
ture dependence of the resistance and from FEM sim- 
ulations that the maximum temperature of the engine 
beam is 370±20 K at I c i c = 3 mA (see Supplementary 
Fig. C). This temperature increase of 17% by DC re- 
sistive heating is relatively small compared to the factor 



5 decrease in T e g as a result of the reduction of Q c g. 
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Supplementary Information 
Piezoresistive heat engine and refrigerator 

In this Supplementary Information the thermodynamic coupling between stress and strain in piezoresistive ma- 
terials carrying a current is discussed in more detail. In Supplementary Discussion A an analytical expression is 
derived for the effective complex Young's modulus Y* s . This expression is used to provide an analytical estimate 
for the feedback coefficient /?, which expresses the effect of a DC current on the effective spring constant k* g of 
a piezoresistive spring. The analytical result is compared to measurements presented in the main manuscript. 
In Supplementary Discussion B the optimal efficiency and power density of silicon heat engines is estimated. In 
Supplementary Figure C a measurement of the resistance of the device as a function of DC current and chuck 
temperature is used to estimate the temperature rise in the engine beam. 



SUPPLEMENTARY DISCUSSION A: Quantitative discussion of the feedback mechanism and coefficient /3 

To evaluate the effect of a DC electric current on the mechanical properties of a piezoresistive solid, consider a solid 
through which a current density Jd c is flowing in the y-dircction. When a small uniaxial vibrational stress <7 ac e lwt 
with frequency uj is present along the y-direction the piezoresistive effect will induce an AC resistivity change p ac : 

Pace 1 " 1 = Pdc ^ace lu)t (Al) 

Where pd c is the unstressed resistivity and 717 is the longitudinal piezoresistive coefficient. It is assumed that 
Pac *C Pdc- This resistivity change modifies the AC resistive heating power density p pr : 

p pr e iut = lz J 2 dc Pa C e luJt (A2) 

As discussed in the main manuscript, jz is a correction factor which is needed if a finite impedance to ground is 
present parallel to the piezoresistive resonator, such that AC currents i ac can also contribute to p pr . From the heat 
equation p = — fc„V 2 T + c p pddT / dt it follows that the heating power causes sinusoidal temperature fluctuations with 
amplitude T ac : 

iait 

T ac e^ = (A3) 

The specific heat capacity is given by c p and pd is the mass density. The effect of thermal heat conductivity kh is 
captured by the factor 7„: 

lh = -k h (V 2 T ac )/T ac (A4) 

The factor jh depends on the resonator's geometry as will be discussed below. The temperature increase generates a 
thermal expansion stress: 

o ac ,tc = a te YT ac (A5) 

Y is Young's modulus and at e is the thermal expansion coefficient. This thermal expansion stress adds to the externally 
applied stress a a c,ext, such that the total stress is given by a ac = a ac .cxt + facte, and the corresponding strain is given 
by £ac = Cac/Y . The stress-strain relation of the piezoresistive solid is therefore the same as that of a solid with an 
effective complex Young's modulus Y* s = <J a c.ext/ £ac- The variables, coupling mechanisms and multiplicative factors 
in equations (|AHIA5P have been schematically shown in figure Al and using these equations the effective modulus can 
be expressed as: 

r; ff H = laazfuE*! „ y (i _ lz PMll3^C \ (A6 ) 

£ac \ Jh + lUpdCp J 

This relation is only valid for small stress and strain, such that the linear approximation is valid and material 
parameters are constant. When all constants in equation (|A6|) are real and positive, equation (|A6|) is identical to the 
Young's modulus of the standard anelastic solid^, however if some of the constants are negative or have a non-zero 
imaginary part, it can lead to properties like mechanical self-amplification and negative creep. This modification of 
the effective Young's modulus can occur in any piezoresistive solid in the presence of electrical current. Besides its 
effect on the dynamics of mechanical structures it might thus also affect the propagation of acoustic waves in solids. 
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elasticity 

strain xY stress + a acm 
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Supplementary Figure A: Schematic of the thermodynamic coupling mechanism via the piezoresistive, resistive heating 

and thermal expansion effects. Relevant variables in the mechanical, electrical and thermal domain are shown. Coupling 

mechanisms are represented by arrows. 

As follows from appendix A of Nowick and Berry^, the effective spring constant of a resonance mode with a complex 
position dependent Young's modulus is given by: 

Ks = k { J v YeUz)dz ) =k{1 ~ (A7) 



Where the integrals run over the volume V of the resonator. Combining this with equation (|A6j) . it is found that the 
feedback coefficient (3 can be expressed as: 

& = Iklz A2 l ? dclTiateY \ = IzPo (A8) 

Where A is the cross-sectional area of the piezoresistive conductor, and 7^ is a factor which accounts for the non- 
uniformity of Y* s . If Y* s is constant over the whole volume, 7& = 1. 

For the resonator under consideration, the thermodynamic feedback effect is large inside the narrow engine beam 
where the strain and current density are concentrated. It is therefore a good approximation to assume that Y* s is 
given by equation (|A6|) inside the volume V^ ng i no of the engine beam and is equal to Y everywhere else in the resonator. 
In this case the factor 7^. is given by: 

L . e 2 ac (z)dz 

iu = Jv y- a ;\; (a 9 ) 

For the lowest frequency in-plane bending resonance mode this fraction is found to be 7^ = 0.11 using a finite element 
simulation. 

To evaluate 7^ the AC heat equation is to be solved: 

khV 2 T ac + p pr = iCpp d uiT ac (A10) 



Since the dimensions of the engine beam are much smaller than the thermal wavelength A/j — y^Sn 2 kh / (c p pd^o) ~ 
26 /im, the heat equation can be simplified by assuming that the AC resistive heating density p pr is zero outside the 
engine beam and that p pr and T ac are independent of position inside the engine beam. The beam can be treated as 
a point source at y = with temperature T ac ^e' lujt . The corresponding solution of the 1-dimensional heat equation 
outside the beam is: 

T(x, t) = Tacfi e- {1+%)2 ^ /Xh+lut (All) 

Since the heat conduction from the 2 ends of the beam needs to be equal to the difference between the generated and 
stored energy in the beam, it follows from equation (|A3[) that: 



- 2fc/iVT ac | y=0 = L(p pr - iCpPdUT ac ) = L"f h Tac (A12) 

Where L=800 nm is the beam length. Substitution of equation (|A1 1|) yields: 

7h = 4tt(1 + i)k h /{L\ h ) (A13) 

The material constants along the [100] direction of n-type silicon are Y — 130 GPa, 717 = — 102 x 10" 11 Pa, a te — 2.6 x 
10 -6 K -i pd =2 329 kg/m 3 , c p = 702 J/kg-K and k h =113 W/(K-m). The cross-sectional area is ,4=280x1500 nm 2 . 
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Substituting these values in equations (|A13|) and (|A8|) and using 7fc=0.11 from equation (|A9|) yields (3q = —132 + 157i 
A -2 . This value is higher than the experimental value, which is mainly attributed to the fact that the actual heat flux 
outside the beam is not 1-dimensional but radial, which results in a larger value of 7^. A finite element simulation of 
the full geometry, including the anisotropy of the silicon crystal results^ in /3q = — 111 + 58i A~ 2 , in good agreement 
with /3q = —123.7 + 64. 6i A -2 as determined from the fits of the measurements in figure 4 of the main manuscript. 



SUPPLEMENTARY DISCUSSION B: Efficiency and power density of silicon heat engines 

The optimal efficiency of any heat engine rjhe,max is limited by the Carnot efficiency r\c — 1 — T m i n /T u 
heat engines there is another efficiency constraint related to the maximum strain before fracture of the solid working 
substance. To illustrate this, consider a silicon beam with cross section A and length L going through an idealized 
Stirling cycle of heating at constant volume, expansion at constant high temperature T max = T min + AT, cooling 
at constant volume and contraction at constant low temperature T m ; n . The beam's spring constant is k = YA/L 
and the static thermal expansion is Axt e = &t e ATL. The heat needed to raise the temperature of the beam is 
Qh = ALc v pdAT , where c v is the heat capacity at constant volume. The isothermal expansion over a distance Ax 
cost external work Wheat, e = ~\k(/S.x — Ax te ) 2 . To keep the temperature stable during expansion, additional heat 
~ Qh(c P — c v )/c v is needed. Then the beam is cooled at constant volume to temperature T m ; n . During isothermal 
contraction, the beam generates external work W CO oi,e = ^fc(Ax) 2 . The total external work done by the engine 
beam is therefore W e = Wheat, e + W CO oi,e ~ kAxAxte = at e YAAxAT, for Ax te <C Ax. The fracture strain limit of 
silicon which is in the ranged of esi.frac ~ 0.01-0.03, as a conservative estimate we use esi.frac — 0.01. The maximum 
expansion Ax from compressive to tensile strain is given by be Ax max = 2egi ! f rac L. Therefore the optimal efficiency 
of a silicon thermal engine is estimated to be: 

rysi, max = W e /Q h « 2 rQte£Si ' frac = 4 x lO" 3 (Bl) 

Pd c p 

Higher efficiencies are possible if the strain is always kept compressive to enable higher Ax, or by using a regenerator 
to reuse the heat lost during cooling. Note that for quasi-static thermal actuators Ax = Axt e and the work done is 
W e = ife(Ax te ) 2 , such that the optimal efficiency is lowered by a factor a te AT/(4£gi,frac)- Even for AT=100 K, a 
quasi-static actuator will have an optimal efficiency that is a factor 150 lower than the optimal dynamic effiency. For 
most quasi-static actuators the efficiency will be much lower because the thermal conduction rate Qt will not only 
cool during the isothermal contraction, but during the whole cycle. 

Now we continue to estimate the maximum efficiency of the piezoresistive heat engine. From figure 5a in the main 
manuscript it is found that the maximum rate of work done by the heat engine is given by Whe = \k < x 2 > ojo^ c Im/3. 
The value of /3 in equation (|A8[) is maximal for n-type silicon (71"; < 0) with 7^ = 1, 7^ = 1. Optimally 7^ is zero, 
however this would mean that there is no thermal conduction and no cooling of the heat engine. As a compromise 
7/j = LopdCp is taken, reducing the efficiency by a factor 2. Thus it is found from equation (|A8[) that the maximum 
value of Yvnfi for the presented piezoresistive heat engine is given by: 

Im /Vma* » - Pd f 2 iateY = 7500 A- 2 (B2) 

The highest measured value of Im (3 is a factor 75 lower, of which a factor 9 can be attributed to -fk = 0.11, 
and the remaining factor ^8 is mainly due to a too high value of 7/j. The heating power in the beam is given by 
Q r = I 2 c pd c L/A. The spring constant is given by k = YA/L and the average amplitude < x 2 >= ^(ssijracL) 2 . Thus 
the optimal efficiency of the piezoresistive heat engine is found to be: 

??pr,Si,max = ~ Z = 3 ' 4 X 10 V B3 J 

A small value value of 7& does not reduce the efficiency, because the reduction in j3 is compensated by the higher 
value of stored energy ifc < x 2 > in the equation for Whe- In other words, the efficiency of the heat engine beam 
is independent of the resonator to which it is coupled. For the heat engine reported in the main manuscript, the 
maximum strain in the engine beam is estimated with finite element simulation to be e& w 0.0015 when operating 
at an amplitude of the center of mass x =100 nm. Moreover the measured Im (3jk is a factor 8 lower than optimal, 
such that the estimated efficiency is ?7he,est = ?7 P r,si,max x (0. 0015/0. 01) 2 /8 = 10~ 6 close to the measured value. The 
difference between the measured efficiency and equation (|B3|) suggest that there is still room for efficiency improvement 
by a factor 500 of the presented heat engine concept by increasing its amplitude and reducing the heat conductivity 
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parameter 7/j. The equations for maximum efficiency (|B3|) and (|B1[) , which are based on linear approximations, 
cannot exceed the Carnot efficiency and should obey rjc > ?7si.max > ?7pr,Si,max- 

Based on the heat equation of a ID beam with fixed temperature at its ends, the DC temperature rise in the engine 
beam is estimated as AT max = -$f- The maximum power density of the piezoresistivc heat engine is given by: 



W-hcmax _ Qr _ ^ Y 2 Q^4,frac k h AT max _ in JW 

Ppr. dens. max — xr — ^?pr.Si.max Tr ~ T n ~ tUU ~ V^^J 

V V p d c p L z m 6 



Where the equation has been evaluated at L=800 nm and AT max = 200 K. Even higher power densities can be 
achieved by reducing the beam length L. When decreasing L the frequency u needs to go up to maintain a sufficient 
efficiency according to the condition 7^ upd,c p and equation (|A13|) . A good heat sink is needed to keep the ends of 
the beam at low temperature. 

Recently much progress has been made on engines operated by electrostatic forces^—. Although these engines can 
have a high efficiency, their work per cycle is limited by the maximum energy^ that can be stored in the electric field 
(.Ee.max «4xl0 5 J/m 3 ). Assuming a 100% efficiency their power density is approximately limited by p e i,dens,max = 
-Ee.max x /■ Comparing this with equation (|B4|) it is estimated that the power density of piezoresistive engines with 
a length L = 800 nm will exceed that of electrostatic engines up to frequencies of 200 MHz. By lowering the length 
L this limit can be pushed to even higher frequencies. Despite their low efficiency, it is thus anticipated that in most 
cases piezoresistive heat engines will outperform electrostatic engines in terms of power density. 



SUPPLEMENTARY FIGURE C: Temperature and current dependence of resistance 



R(Q) 




l«>A 2 ) 



Supplementary Figure C: Measurement of the resistance of the resonator Rd c as a function of Idc and chuck temperature 
Tchuck- The difference in resistance between two equiresistance lines is 1.4 Q. For a uniformly heated device, the temperature in- 
crease AT is proportional to the DC power Rdcldc an d the slope of the equiresistance lines in the figure is expected to be constant. 
In the measured temperature range, the slope of the equiresistance lines is between AT//J c =0-8 K/(mA 2 ). A FEM simulation 
of the structure with a constant (temperature independent) resistivity of ptj c =lCP 4 f2m predicts AT//J C =3.75 K/(mA 2 ) for the 
maximum engine beam temperature^. The observed variation in the slope is due to the non-uniform temperature distribution in 
the structure in combination with the strong temperature dependence of TCR, which is near a minimum and is even observed to 
be slightly negative at room temperature—. As a consequence of this, the total resistance of the structure is almost independent 
of Idc around T c h uc k=30 o C, where the increase of resistance of the warmer parts compensates the decrease of resistance of the 
colder parts of the structure. From these data it is estimated that in practice AT//J C =6±2 K/(mA 2 ), such that the maximum 
temperature rise in the beam is estimated to be AT=56±18 K at 3 mA at maximum cooling power and AT=7±2 K when 
operating as a heat engine at 1.05 mA. 
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